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ABSTRACT 



Aims. We report a simple method to generate potential/surface density pairs in flat axially symmetric finite size disks. 
Methods. Potential/surface density pairs consist of a "homogeneous" pair (a closed form expression) corresponding to a uniform disk, 
and a "residual" pair. This residual component is converted into an infinite series of integrals over the radial extent of the disk. For a 
certain class of surface density distributions (like power laws of the radius), this series is fully analytical. 

Results. The extraction of the homogeneous pair is equivalent to a convergence acceleration technique, in a matematical sense. In the 
case of power law distributions (i.e. surface densities of the form T.(K) oc R s ), the convergence rate of the residual series is shown to be 
cubic inside the source. As a consequence, very accurate potential values are obtained by low order truncation of the series. At zero 
order, relative errors on potential values do not exceed a few percent typically, and scale with the order N of truncation as l/N 3 . This 
method is superior to the classical multipole expansion whose very slow convergence is often critical for most practical applications. 

Key words. Gravitation — Methods : analytical — Accretion, accretion disks 



1. Introduction 

The construction of potential/density pairs i s a common concern 
in the context of galactic dynamics (e.g. IBinney & Tremaine 
Il987l) . More generally, the accurate determination of the grav- 
itational potential corresponding to a given mass distribution 
is a critical step in astrophysical disk models and numeri- 
cal simulations. Unfortunately, potential/density pairs are not 
closed-form expressions, except in a few cases only dMestell 
Il963t lde Zeeuw & Pfenniger|[T988l: IBinney & Tremainelll987l) . 
and most known pairs involve infinite ma t ter distribution (e.g. 
Ouar3fl99H ICuddefordlfl99H llarnl[T996t IBinney & Trem aine 



1987b . whereas astrophysical disks are finite (in size and mass). 



For flat (i.e. zero thickness) and finite disks of interest here, 
the potential is mainly accessible through the double integration 
of th e surface density weighted by the Green function l/\r - r'\ 
(see IBinney & Tremaing|l987l for a review). This integral ap- 
proach is very tricky due to the inevitable occurrence of singu- 
larities everywhere inside matter, i.e. when r — > r' '. The dif- 
ficulty is c ircumventet w hen the Green function is converted 
into series (Kellogg 1929). However, these are infinite and alter- 
nate series which do not converge insid e sources, alth ough their 
integration over the distribution does dDurand| [T953). The effi- 
ciency of the series representation is therefore strongly mislead- 
ing in practice because series truncations do not lead to a numer- 
ically stable description of the potential inside disks dClementl 
119741) . The errors in potential values can be as large as a few 
percent on average, even when the order of ttuncation attain 
several tens dStone & No rman 1992). Singularities are no more 
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present when the potential is expressed in terms of integrals of 
Bessel functions. But the oscillatory behavior of Bessel func- 
tions as well as their infinite definition domain pose technical 
difficulties, even for infinitely extended disks. Potential/density 
pairs can also be constructed by s uperposition of homeoid s con- 
veniently shrunk along one axis (Binn ev & Trem aine 1987|), or 
by param etric expansion and or dering of the associated Poisson 
equation dCiotti & Bertinll2005l) . 

iHure & Pierensl d2005l) have reported a numerical method 
to compute the gravitational potential in flat axially symmet- 
ric disks with great accuracy by a correct treatment of the sin- 
gularity. When the surface density profile is conveniently split 
into two components (namely, a homogeneous component plus 
a residual one), the contribution of the singularity to the potential 
can be accounted for exactly through a closed-form expression. 
In this paper, we explore this "splitting method" from a purely 
theoretical point of view in search for fully analytical solutions 
(i.e. potential/surface density pairs). In particular, it is shown that 
the residual component is an infinite series of radial integrals 
which are known analytically for a certain class of surface den- 
sity distributions. This is the case of p ower law pro files which 
are of great astrophysical interest (e.g. |Pringld [l981). The con- 
vergence rate of the residual component is then cubic. This is a 
major issue, when compared with the classical multipole expan- 
sion whose converge speed is prohibitively low inside sources. 

This paper is organized as follows. In Sect. [2] we briefly 
recall the density splitting method for one dimensional disks 
(Hur e & Pierensl 120051) . We expand the residual component as 
a infinite series of the integrals over the modulus of the com- 
plete elliptic integral of the first kind. We then show that the 
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construction of potential/surface density pairs is possible for a 
certain class of surface density distribution. The impact of the 
series truncation on potential values is discussed in Sect. [3] The 
case of surface densities scaling as a power of the radius is in- 
vestigated in detail in Sect. [4] (a Fortran 90 code is available on 
request). For this kind of distribution, the convergence rate is 
cubic. This is established in Sect. [5] The paper ends with a sum- 
mary and concluding remarks. 




Fig. 1. A finite axially symmetric disk. 



2. Exact potential/surface density pairs in flat 
axisymmetric disks 

2. 1 . Outline of the density splitting method 

The gravitational potential *F in the plane of a fiat axially sym- 
metri c disk as pictu red in Fig.[T]is given by the exact expression 
(e.g. lDurandlfl95l : 



= -2G f 



2(a) A / —mK(m)da, 



where 



K(m) = 



n/2 



d(p 



■ m 2 sin 2 < 



is the complete elliptic integral of the first kind, 



a+R 



(1) 



(2) 



(3) 



is the modulus of K (with < m < I), a and R are polar cylin- 
drical radii (a refering to the matter distribution), a; n > is the 
radius of the inner edge, a out > a m is the radius of the outer 
edge, 2(a) is the surface density, and G is the Gravitation con- 
stant. The function K is plotted versus m in Fig. |2j it is char- 
acterized by a dive rgence for m — > 1 . This is a logarithmic di- 
vergence (see e.g.lKellogd | 1929tlMeste]| [T963), since we have 
(Gradshtevn & Rvzhik 1965): 



lim K(m) = In 



[l +0(1 -m 2 )]. 



(4) 



Despite this singularity which occurs everywhere inside the disk 
(i.e. for a m < R < a out ), the potential is generally finite every- 
where. 

In order to determine *F inside the disk, we set: 



u = — < 1 
R 



and 



v=*<l, 
a 



(5) 



and so, u locates matter of the inner disk [a m ,K\, whereas v 
is for the outer part \R, ow l. Using the transformation (e.g. 
iGradshtevn & Rvzhikll 19651) : 



K|^| = (l + .v)K(.v). 
1 + x 



(6) 
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Fig. 2. Complete elliptic integrals of the first and second kinds 
versus the modulus m. 



V(R) 



-4GR 



(7) 



where < x < 1, Eq.(Q~|l reads: 

f l.(u)K.(u)udu — f S(v) — ^-dv 
J«i„ Ji v 1 

where w m = a m /R < 1 and v out = R/a out < 1 correspond to 
the position of the disk edges. As quoted, S must now be re- 
garded as a function of u in the first integral, and as a function of 
v = 1/u in the second one. Ac cording to the "density splitting" 
method dHure & Pierensl l2005h . we divide 2 into two compo- 
nents, namely: 

2(a) = E +SUa,K), (8) 

where 2o = 2(/?) is the value of the surface density at the radius 
of the singularity (where the moduli m, u and v are unity), and <52 
is the "residual" surface density profile. The potential then takes 
the form: 

TO = V (R) + ff¥(R), (9) 

where 



V (R) = -4G/?2 



J" K.(u)udu — J" 



K(v) 



dv 



(10) 



is the potential due to a radially homogeneous disk with constant 
surface density 2o, and 



ff¥(R) = -4GR 



f 



6'L(u)K(u)udu 



I 



K(v) 
5Z(v)-Y-dv 



(11) 



is the "residual" potential associated with 52. Note that 2o varies 
from place to place in the disk , as well as <52. As demonstrated 
elsewhere (see also Hure 2005), this technique of density split- 
ting can be numerically very accurate, for two reasons. First, *Fo 
is known in a closed form: 



W) = -4Gfl2 



E(v out ) t 2 

- E(M in ) + u- K(u in ) 



^out 



(12) 



and is finite everywhere in the disk mid-plane. In Eq.(fT2l. 
u 'ia ~ -^1 - u 2 n is the complementary modulus and E is the 
complete elliptic integral of the second kind (see Fig.|2|: 



E(m) 



m 2 sin (p. 



(13) 
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Second, each integrand in Eq.dTTb now vanishes at the singu- 
larity (this is the role of the splitting), making the two integrals 
easily computable with standard quadrature rules. 

2.2. The residual potential from series representation 

In order to determine analytical expressions for Y on the 
basis of the splitting method, we consider the expansion of 
the elliptic integral of the first kind over its modulus (e.g. 
iGradshtevn & Rvzhikll 19651) . namely: 



K(m)=^£y„m 2 ", 



To = 1, 



(14) 



In = 



(2n- 1)! 



2"n! 



= 7n-l 



2n-l 
2n 



n>\. 



When this infinite series is inserted into Eq.dTTb. the residual 
potential becomes an infinite series. Let 



n=0 



be the sum of N + I terms xo, x\, . . . , x^. The residual potential 
is then given by: 

ff¥ = S 00 (6V„), (16) 



with for any n > 
<5¥„ = -2nGRy n 



f 6T.(u)u 2n+l du 



- ^ °" 5I.(v)v 2 "- 2 dv 



(17) 



Note that this expression for 6^ is still exact, provided all the 
terms of the series are included. It follows that poten- 
tial/surface density pairs (Sfo + So + SI.) can easily be con- 
structed for flat disks with finite size (or not, depending on u m 
and v out ). These pairs are analytical if the residual pair (5*T, 51.) 
is analytical, that is, if the product of E(a) by a power law of a 
is analytically integrable. We immediately see that many simple 
functions £(a) are concerned. In particular, if X is a power law of 
the radius a — as often met in astrophysical disk models — , then 
d^V,, is a power law of the radius too (more precisely, a mixture 
of four power laws; see Sect. 0). Obviously, pairs determined 
through Eqs.dTSb and ( fT7b should generally remain in the form 
of an infinite series with no equivalent closed form expression. 
This is not a big problem, first because working with infinite ex- 
pansions is common in potential theory, and second because the 
residual potential is a rapidly converging series (see Sect. [5]). 



3. Approximate potential from truncated series 

3.1. General remark 

The estimate of should generate no significant error since 
special functions E and K can be determined within the com- 
puter precision from a numerical library. In contrast, the com- 
putation of d^V is not errorless because, in practice, any series 
is necessarily truncated at a certain order N. Does the trunca- 
tion produce large errors in potential values? In order to answer 
this question, it is important to stress again that the occurrence 
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exact (no approx., 
one term expansion of K 
two term expansion of K 



0.2 0.4 0.6 0.8 1 

modulus u and v 

Fig. 3. Typical shape of the integrands appearing in Eq. dTTb com- 
puted from different approximations for the function K. In this 
example, the disk has inner edge a m - 0.1 and outer edge 
flout = 1, the surface density obeys a power law of the radius with 

>(flin"t" flout)- 



( 15 - ) E(a) oc a ~ 3/2 (that is s = -1.5), and the radius is R 



of the condition SE — enables the neutralization of the loga- 
rithmic singularity, making the remaining integrands fully regu- 
lar. As a consequence, the value of the residual potential 6^ is 
mainly determined by the shape of the integrands far from the 
singularity, that is close to the edges. In other words, neither 
the precise behavior of the complete elliptic integral K nor the 
shape of the residual surface density profile 61. around the sin- 
gularity are critical to estimate r5 v F with a certain accuracy. We 
therefore expect that 6*¥ can be accurately determined if K is re- 
placed by any approximation which preserves at best the values 
of this function around the bounds u m and v out rather than at the 
radius of the singularity (where m = u = v = 1). From this point 
of view, the series defined by Eq.dTTb is fully appropriate as it 
converges very rapidly for moduli significantly less than unity. 
Figure [3] shows a typical example where S is a power law of the 
radius. We clearly see that the two integrands 6"L(u)uK(u) and 
(5E(v)K(v)/v 2 appearing in Eq.dTTb are quite insensitive to the 
approximation adopted for K (one-term, two-term expansions 
or exact). This conclusion is independent of the Z-profile and of 
R as well. 

3.2. The first terms 

For a truncation of the series at order N, the potential is given by 
the approximate value: 



¥ approx . = ¥o+ > WJ'„)*Y > 
and we have asymptotically: 



¥ + lim [S N (ff¥ n )] = ¥. 



(18) 



(19) 



With only one term in Eq.lTBb. we have 5o(5T'„) = 5*¥o with: 



-2nGR 



p 1 r ,w dv 

I 61.(u)udu - I (5£(v) — 
J«,„ Ji v 2 



(20) 



and so the approximate potential is Yapp,™. = To + £*IV For a 
two-term expansion of the complete elliptic integral of the first 



4 



Jean-Marc Hure, Didier Pelat and Arnaud Pierens: Potential/surface density pairs in flat disks 



kind, the residual potential is S i (£*P„) = 6^0 + S*¥i , where 6^0 
is given by Eq.d2"0"l), and: 



<5T] = — GR 
2 



r 1 •? r Vmt 

I 6L(u)u du - I £E(v)aV 



(21) 



The approximate potential is then *P a 



and so on. Each new term added into the expansion brings a 
new component <S*F„ defined by Eq.dTTb. We then expect that the 
more terms in the expansion, the more accurate the potential. If 
the series converges rapidly (this is the case of power law dis- 
tributions; see below), the absolute error made on *P is roughly 
l<^V,|. 

4. The case of power-law surface density profiles 

Let us consider the case where 2 is a power-law of the radius, 
namely: 

a 



S = cr 

\ flout 



(22) 



where <x is the surface density at the outer edge of the 
disk. Power-law s urface density solutions are found in many 
disk models (e.g. Pringle 1 198 lh and match observations (e.g. 
iGuilloteau & Dutrevll998 ). with s < in general. We then have: 



and so: 



61. : 



So (v~ 



1) for u < 1, 
-1) for v<l. 



(23) 



(24) 



4.1. Order zero 

For a one-term expansion of the complete elliptic integral of the 
first kind, the residual potential is given by: 

^ ^ f* v ° ut dy 

{\-u s )udu+ | (v~' 5 - 1) - 
Ji V 

-1), this is (see AppendixlAlfor these two cases): 
i s+2 - 1 1 - u 2 

m , m 



6% } = InGR^) 



f 



(25) 



(5^0 = 2nGRZ l 



s + 2 



l-v; 



1 



s + 1 



(26) 



Figure |4] displays the approximate potential ^pp,™. = 
6*¥o for disks with inner edge «;„ = 0.1, outer edge a out = 1 and 
power law exponents s = — 1, — 1.5 and -2.5. References poten- 
tial values *F re f (R) have been determined withi n the computer 
precis ion from the numerical splitting method dHure & Pierensl 
120051) and compared to Tapprox.- The decimal logarithm of the 
relative error e = |1 - Tapprox./^ief.l is a l so shown in Fig. [4] We 
see that the potential is reproduced with an accuracy better than 
~ 2% (on average) in the entire disk with only one term in the 
expansion of the K-function. This is already remarkable. The ac- 
curacy is relatively uniform, and it is slightly better at the outer 
edge. This point is particularly interesting to model young stel- 
lar objects or Active Galactic Nuclei. In these systems, the inner 
disk is dynamically under the control of a massive central object 
(a proto-star or a black hole), and so there is no need for very 
accurate potential values. On the other ha nd, in the ou ter regions 
where disks can be self-gravitating (e.g. Hure 2000), a reliable 
description of the potential is required. 
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Fig. 4. Top: potential (lines) inside the disk for three different 
surface density exponents s when computed with a one-term ex- 
panded elliptic integral (i.e. S'V = S^o). Reference values T^ef. 
(symbols) are obtained numerically within the computer preci- 
sion. Bottom decimal logarithm of the relative error (the average 
of e is also indicated). 

4.2. First order 

The accuracy of potential values can be improved by adding a 
second term. For = 1, we find: 



c^i = -GRZ Q 



f 



(1 - u s )u du + 



l)dv 



(27) 



If s + {-4, 1} (see Appendix iBl for these two cases), the ex- 
pression is: 



5*¥i = 2nGRI» 



,s+4 



4(s + 4) 



1 1 -u 4 
+ 



16 

1 - v-" 1 

out 

4(^-1) 



(28) 



Figure|5]displays the results obtained with a two-term expan- 
sion of K, i.e. ^Tapprox = ^o + + under the same condi- 
tions as for Fig. [4] We see that the relative accuracy of mid-plane 
potential values (of the order of 0.3% on average) is now better 
by a factor of ~ 5. Such a level of accuracy is probably suffi- 
cient for many disk models and applications. However, since the 
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Fig. 5. Same legend as for Fig. [4] but the residual potential 5*¥ 
is determined from a two-term expanion of the complete elliptic 
integral of the first kind (i.e. {?¥ = S^q + 8¥i). 



process of adding more and more terms represents no technical 
difficulty and converges rapidly (see Sect. |5j, it can be repeated 
until a given accuracy is reached. In particular, for s + 2n + 2 ± 
and s — In + 1 + 0, we have (see the Appendix ICl for these two 
cases) 



8V n = 2nGLoRy n 



1 — u 



2(n+l) 



1 — U 



j+2(n+l) 



2(n+l) s + 2(n+l) 

v 2n-s-\ _ J 



v 2 "; 1 - 1 

*out 



2n — s — 1 



2n- 1 



(29) 



5. Convergence rate of the series 5*F = 5 in 
the case of power law distributions 

5.1. A cubic convergence ? 

Figure [6] shows the rise of the accuracy of potential values ¥(/?) 
as the order N of the truncation increases for the same disk pa- 
rameters as already considered and a few exponents s. We find 
that the relative error e = 1 1 — v P a pproit,/ l Pref,l fi rst decreases slowly 
as N increases. But above a certain rank «o ~ 8 - 10, this error 
approximately scales as 1 /TV 3 , before saturation near the com- 
puter precision, for N > 10 4 terms. We have checked that this 
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Fig. 6. Decimal logarithm of the relative error on potential values 
*F = 1*0 + versus the order N of the truncation in Eq.(fT9b for 
five power law exponents s = -3.5, -3, -2.5, —1.5 and -0.5. 
The disk has inner edge a m = 0.1, outer edge a out = 1 and the 
potential is computed at R = l(a\ n + a out ). 




10 100 
order n 



10 100 
order n 



Fig. 7. Left: y n versus n. Right and r„ (for s = -1.5) versus n. 



conclusion holds for a wide variety of configurations and radii 
R. Convergence appears to be cubic. With this approach, we can 
easily select the lowest order corresponding to a given level of 
accuracy. From Fig. [6] we have roughly: 



loge « -31ogN-2±0.5, 



(30) 



for N > n , or N ~ (100e) _1/3 within about 0.17%. 

There are two remarkable cases. For 5 = 0, the potential is 
fully analytical: *P = (i.e. = 0). For s = -3, the relative 
accuracy rises with N much more drastically than for others ex- 
ponents and reaches the maximum allowed by the computer for 
N > «o (see below). 



5.2. Proof 

The expansion of K(ra) through Eq. ( fT4b converges slowly as m 
approaches unity, and even diverges for m = 1, although y„ de- 
creases rapidly as n increases (for large n, y„ ~ I In). This is 
shown in Figs.|7]and[8] Note that the divergence of K as m — > 1 
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Fig. 8. Left: series S /v(y«) versus the truncation order N which 
diverges as N — > oo. Right: series 7tSai(£„)/4 versus the trun- 
cation order N which converges towards Catalan's constant (see 
also Fig. [7J. 



is the main reason why the determination of the potential in- 
side sources by direct integration is generally not the appropri- 
ate way. What about the convergence of 5*¥ ? Since 5*¥„ comes 
from the integration of K over its modulus, S aK^m) is made of 
terms of the form ^j-j = £ n which decrease faster with n than 
terms y„ do. The series d^Y is then expected to converge (see the 
Appendix iDb. Figure [7] displays the terms £„ versus n, and the 
series S n(£„) is plotted versus N in Fig. [8] 

When the series r5*P is truncated at order N, the relative error 
e(N) on the potential value is: 



A approx. A 



From Eqs.(fT8b and il9i . this expression can be re-written in the 
following form: 



e(N) 



(32) 



The quantities u m and v out being always smaller than unity, 
their contribution in Eq.d29l) — an exponential drop — is ex- 
pected to vanish rapidly as n exceeds a certain rank. This tran- 
sition occurs at Af ~ no; it is clearly visible in Fig. [6] where 
no ~ 8. Thus, for N > no, the leading terms in d 1 ? form a series 
SivCr,,), with: 



Jn 



1 



1 



1 



2(«+l) 2(«+l) + s In - 1 
1 



(33) 



2/1—1 — 5 

which can easily be rearranged: 

s(s + 3)(4n+ 1) 



T " Jn {In + 2)(2n + 2 + s){2n - l)(2n - 1 - s) ' 
For large values of n, we have: 



(34) 



(35) 



and so, r„ behaves like n~ 4 , asymptotically. The coefficient r„ 
is displayed versus n for s = -1.5 in Fig. [7] Then, as soon as 
N » «o, the relative error is 



e(A0 



2nGI.oR 



I N 



(36) 



Using Eq.d35ll, we have: 

2nGT.oR 



T„dn 



27tGZoRs(s + 3) 



Jn 



(37) 



dn 



¥ Joo 4n 3 

Since y„ scales as 1 jn for large n (see also Fig.©), we have: 

e(N)oc — , (38) 

as found "experimentally". The convergence rate of this series 
is therefore cubic. This result holds irrespective of the radius R 
(except at the disk edge£3) and whatever s, except in the case of 
the homogeneous disk (i.e. for s — 0) and for s = -3 as well. In 
this latter case, the convergence rate is not limited by the terms 
t„, but by the series asssociated with w m and v out . For s = -3, we 
have: 



8¥ n = 2KGY.oRjn 



2(«+l) 



U 



2(n+l) 



„2n- 



2(n+ 1) 



2n- 1 



(39) 



and so the convergence of &¥ is exponential (see Fig. [6j, except 
at the inner and outer edges (see Appendix|E]i. 

All numerical experiments have succesfully illustrated these 
results. This is a major issue, in particular because classical mul- 
tipole expansions converge very slowly inside sources. 



(31) 6. Summary and concluding remarks 



In this paper, we have reported a method to determine exact and 
approximate potential/surface density pairs in flat finite size ax- 
ially symmetrical disks, in the continuity of the numerical split- 
ting method described in Pierens & Hure (2005). At a given ra- 
dius R, each pair is the sum of two components: 

- a "homogeneous" pair (T^, ^o) where Eo = £(/?) is the local 
value of the surface density and To is given by Eq. (fT2l . 

- a "residual" pair (<5 V F, SI.), where 5*? is an infinite series (see 
Eqs.© andO) and 61(a) = 1(a) - £ . 

In general, the infinite series 6^ cannot be put into a closed 
form expression. However, in contrast to the multipole expan- 
sions commonly met in potential theory, this series has a cubic 
convergence rate inside the source. This major result is due to 
the splitting method which is equivalent to a convergence accel- 
eration technique (in the mathematical sense). It means that very 
good approximations for the potential can be obtained by low 
order truncations of the series. Pairs can be constructed only for 
a certain class of the surface density profile 2(a): any function 
2(a) whose product by a power law of the radius a is analytically 
integrable does lead to a potential/surface density pair. This is the 
case of power laws. It would be interesting to list all functions 
satisfying this criterion. 

We have fully considered power law surface densities. Such 
profiles are particularly attractive for two main reasons: i) power 
laws underly many (if not most) astrophysical disk models, and 
ii) power laws may serve as a basis set to construct more com- 
plicated surface density profiles. This new approach along with 
associated computing toolfl should enable a much better repre- 
sentation of the gravitational potentials and forces in flat disks. 



1 A similar calculus shows that convergence is only quadratic. 

2 An easy-to-use Fortran 90 package nammed APPLawD (an acronym 
for "Accurate Potential for Power Law Disks") is available upon request 
to the authors. APPLawD can be downloaded at the following address: 
www . obs . u-bordeauxl . f r/radio/ JMHure/ intro2applawd . html 
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Appendix D: Convergence of S^n) 

From (e.g. Gradshteyn & Ryzhik 1965), we see that 
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Appendix A: Special cases for dY 

For s = — 1 in Eq.(l22b. then Eq.(l20ji reads 



m 2n dm, 



(D.l) 



(5% = 2nGY. Q R • 



For s — -2 in Eq. (l22l . then Eq.(l20l> reads 



1 


[» 2 


l 








u 


+ 


In v + - 






2 


"in 


V 


1 J 



6%) = 2ttGI,oR- 



- In u 



v + 



(A.l) 



(A.2) 



Appendix B: Special cases for 6Wi 

For s = -4 in Eq.d22b. then Eq.d2"n> reads 



SW t = -GXqR < 



if 

In u 

4 



For s — + 1 in Eq.(l22l. then Eq.d2Tb reads 



u 4 u 5 



+ [\nv-vt 



(B.l) 



(B.2) 



4 5 

Appendix C: Special cases for 8*¥ n 

For s + In + 2 = in Eq.(|22]i, then s - In + 1 + and Eq.dTTTi 
reads 



6¥ n = 2nGLaRy„ 



.201+1) 



2(71 + 1) 



In u 



,,2n-s-l 



(CI) 



„2n-l 



2n — s — 1 2n - 1 

For s + 2n + 2 + and s - 2n + 1 = in Eq.<l2"2l). then Eq.(fT7b 
reads 



W„ = 2nGY. Q Ry n < 



it 



2(n+l) 



It 



s+2(n+l) 



1 



2(n+l) 5 + 2(n+l) 



(C2) 



,,2n-l 



In v — 



In- 1 



n=0 

= 2_i I Jnm 2 "dm, 

^ u n=0 

2 f 1 
= — I ~K(m)dm, 
x Jo 
4C 



where C is Catalan's constant. 

Appendix E: Convergence rate for s = -3 at the disk 
edges 

At the disk outer edge where v out = 1, we have: 
I — a.20>+1) i — A 2 "- 1 



W n = 2nG^ Q Ry n 



2(71 +1) 2ti-1 



(E.l) 



where A = ai n /a ut is the disk axis ratio. Since A < 1, the leading 
terms in 5^¥ then form a series S n(t'„), with: 



Yn 



1 



1 



-Yn 



2(7i+l) 2«-l 
3 

2(n+ 1)(2«- 1)' 



(E.2) 



For large n, we have t' h ~ 1/n 3 , and so e(N) oc 1/N 2 . The con- 
vergence is therefore quadratic (as at the disk inner edge). 



